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Wind structure of late B supergiants. 
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ABSTRACT 



Aims. We provide a quantitative analysis of time-variable phenomena in the photospheric, near-star, and outflow regions of the late-B supergiant 
(SG) HD 199 478. This study aims to provide new perspectives on the nature of outflows in late-B SGs and on the influence of large-scale 
structures rooted at the stellar surface. 

Methods. The analysis is based primarily on optical spectroscopic datasets secured between 1999 and 2000 from the Bulgarian NAO, Tartu, 
and Ritter Observatories. The acquired time-series samples a wide range of weak metal lines, He i absorption, and both emission and absorption 
^ ■ signatures in Ha. Non-LTE line synthesis modelling is conducted using FASTWIND for a strategic set of late-B SGs to constrain and compare 
k> their fundamental parameters within the context of extreme behaviour in the Ha lines. 

, Results. The temporal behaviour of HD 199478 is characterised by three key empirical properties: (i) systematic central velocity shifts in the 
' photospheric absorption lines, including C n and He i, over a characteristic time-scale of ~ 20 days; (ii) extremely strong, variable Ho- emission 
with no clear modulation signal, and (iii) the occurrence in 2000 of a (rare) high- velocity absorption (HVA) event in Ha, which evolved over ~ 
60 days, showing the clear signature of mass infall and outflows. In these properties HD 199478 resembles few other late-B SGs with peculiar 
emission and HVAs in Ha (HD 91 619, HD 34085, HD 96919). Different possibilities accounting for the phenomenon observed are indicated 
and briefly discussed. 

Conclusions. At the cooler temperature edge of B SGs, there are objects whose wind properties, as traced by Ha, are inconsistent with the 
predictions of the smooth, spherically symmetric wind approximation. This discordance is still not fully understood and may highlight the role 
of a non-spherical, disk-like, geometry, which may result from magnetically-driven equatorial compression of the gas. Ordered dipole magnetic 
fields may also lead to confined plasma held above the stellar surface, which ultimately gives rise to transient HVA events. 

Key words, stars: early-type - stars: SGs - stars: fundamental parameters - stars: winds, outflows - stars: magnetic fields - stars: individual: 
HD 199478 



1. Introduction 

The key limiting assumptions incorporated within current hot 
star model atmospheres include a globally stationary and spher- 
ically symmetric stellar wind with a smooth density stratifica- 
tion. Although these models are generally quite successful in 
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describing the overall wind properties, there are numerous ob- 
servational and theoretical studies which indicate that hot star 
winds are certainly not smooth and stationary. Most (if not all) 
of the time-dependent constraints refer, however, to O-stars and 
early B supergiants (SGs), while mid- and late-B candidates are 
currently under-represented in the sample of stars investigated 
to date. 
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Vink et al. (l2000h have shown that for mid and late-B 
SGs there is a discrepancy between theoretical predictions 
and Ho- mass-loss rates, derived by means of unblanketed 
model analysis. This finding was confirmed by recent inves- 
tigations using line-blanketed model atmos pheres (see, e.g., 
Crowtheret all 120061: iMarkova & Pulsl 12008b . The reason for 



the discrepancies is not clear yet but wind structure and 
variability might in principle cause them. Indeed, observa- 
tions indicate that while winds in late -B SGs are significantl y 
weaker than those in O SGs (e.g. Ma rkova & Pulsl [2008), 
there is no currently established reason to believe that weaker 
winds might be less structured/variab le than stronger ones (e.g. 
IMarkova et al.ll2005llPuTs et al.ll2006h . 

The first extended spectroscopic monitoring campaigns of 
line-profile variabilit y (Ipv) in late-B SGs were performed by 
iKaufer et al.1 dl996alfbh . who showed that stellar winds at the 
cooler temperature edge of the B-stars domain can also be 
highly variable. Interestingly, in all 3 cases studied by these au- 
thors, the variability patterns (as traced by Ha) were quite sim- 
ilar consisting of (i) blue- and red-shifted emission with V/R 
variations similar to those in Be-stars, and (ii) the sudden ap- 
pearance of deep and highly blue-shifted absorptions (HVAs). 
Though the kinematic properties of the HVAs in Ha were 
found to be completely different from those of DACs (Discrete 
Absorption Components) in the UV spectra of O and early-B 
stars (e.g., HVAs do not propagate outwards, but instead extend 
to zero velocity and even indicate mass infall), similar scenar- 
ios consisting of large-scale wind structures rooted in the pho- 
tosphere were suggested to interpret their appearance and de- 
velopment in time. 

The present paper is focused on a multi-line investigation 
of HD 199478 (HR 8020), a B8 Iae star whose stellar and 
wind properties have been recently determined by means of 
NLT E model atmosphere ana lysis of strategic lines in the op- 
tical (IMarkova & Pulsl 120081) . The first extensive monitoring 
campaign of Ipv in the optical spectrum of this star revealed 
Ha variability similar to those described by Kaufer et al. (i.e., 
peculiar wind emission with V/R variations similar to those in 
Be stars) with one exception, i. e. no indications for any H VAs 
were found during that survey (IMarkova & ValchewOOOh . In 
addition to wind variability significant absorption Ipv was also 
established, which raised suggestions of a link stellar pulsa- 
tions. To check the pulsational hypothesis however long-term 
photometric observations were required. 

Motivated by the intriguing time-variable properties re- 
ported above, we organised and conducted new paral- 
lel spectroscopic and photometric monitoring campaigns of 
HD 199 478 during 1999 and 2000. Details of the spectroscopic 
analysis are presented here and th e results of the pho tometric 
survey were recently published by ( Per cy et al. Il2008l) . 



2. Observations and data reduction 

2.1. Spectral data 

Spectroscopic data consisting of 65 spectra centered on Ha and 
46 on He 1 /15876 were predominantly collected at the National 
Astronomical Observatory (NAO), Bulgaria while individual 



Ha observations were also secured at the Tartu Observatory 
(TO), Estonia and at Ritter Observatory (RO), USA. 

The total time coverage of these data is from January 1999 
to December 2000 with large data gaps in the summer and 
the winter each year. The time sampling was typically 3 to 
6 spectra per month with a time-interval between successful 
exposures of 1 to 2 days, except for the fall of 2000 when 
HD 199478 was monitored more intensively. The time distri- 
bution of the data and spectral regions observed are given in 
Table [2] while specific information about the equipment, reduc- 
tion strategy and methods used at each observatory is outlined 
below and summarised in Table [1] 

National Astronomical Observatory, Smolyan, Bulgaria - 
A total of 95 high-quality spectra (49 in Ha and 46 in 
He 1 /15876) were obtained (observers T. Valchev and H. 
Markov) in the coude focus of the NAO 2m telescope of the 
Bulgarian Academy of Sciences during the period March, 1999 
- December, 2000. Most of the observations (91) were carried 
out using a BL632/14.7 grooves mnT 1 grating in first order 
together with a PHOTOMETRICS CCD (1024 x 1024, 24p) 
as a detector Q. This configuration produces spectra with a re- 
ciprocal dispersion of ~0.2 A pixel -1 and an effective resolu- 
tion of ~ 2.0 pixels (-0.44 A R=15 000) over a wavelength 
range of ~204 A . In addition, four more Ha spectra with a 
reciprocal dispersion of ~0.1 A pixel 1 and an effective reso- 
lution of ~ 2.0 pixels (~0.2 A R=30000) were obtained using 
a BL632/22.3 grooves mmr 1 grating in second order and the 
same detector. The S/N ratios of the NAO spectra range be- 
tween 250 and 450. 

We followed a standard procedure for data reduction (de- 
veloped in IDL) including: bias subtraction, flat-fielding, cos- 
mic ray removal, wavelength calibration, correction for helio- 
centric radial velocity (V r = -12 km s _1 ), water vapour line re- 
moval and rebinning to a step of 0.2 A per pixel. More infor- 
mation about the reduction procedure can be fo und elsewhere 
(|Markova & ValchevlEoOol: IMarkova et alJl2004 . 



Tartu Observatory, Toravere, Estonia - At TO (observer I. 
Kolka) 9 Ha spectra were obtained using a 1.5-m reflector 
equipped with a Cassegrain spectrograph and an Orbis-1 (LN2- 
cooled) CCD camera (512x512 pixels). About 130 A is cov- 
ered by one exposure with a reciprocal dispersion of 0.25 A/pix 
resulting in a resolution of about 16 000 at Ha . 

The observations were reduced in a uniform way using 
MIDAS. Due to the very low dark current, the mean back- 
ground subtracted from the raw frames is a sum of bias plus 
dark and real sky. Flat-fielding was not performed due to the 
reasonably flat response of the CCD and to the empirical re- 
sult that the spectrum summed over 4-6 CCD rows has almost 
no distortion from pixel to pixel sensitivity differences. Thus, 
the photon noise and the read-out noise are the main sources 
of errors. The telluric water vapour lines were removed by di- 
viding individual spectra with a scaled model telluric spectrum. 



1 This detector was characterised by a rms read-out noise 3.3 elec- 
trons per pixel (2.7 ADU with 1.21 electrons per ADU) 
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Observatory Telescope Instrument R-X/8X window [in A] region N s 



National Astron. Obs. 


NAO 


2.0m 


coude 


15 000 


204 


Ha 


45 


National Astron. Obs. 


NAO 


2.0m 


coude 


30000 


100 


Ha 


4 


National Astron. Obs. 


NAO 


2.0m 


coude 


15 000 


204 


He i ,15876 


46 


Tartu Observatory 


TO 


1.5m 


long-slit 


16000 


127 


Ha 


9 


Ritter Observatory 


RO 


1.1m 


fiber-fed 


26000 


70 


Ha 


7 



Table 2. Summary of the spectral data sets. 



Region 


Observational dates 


HJD 245 1200+ 


N 

A y spec 


S/N 


Ha + CII 


1998 Dec. ,30 - May, 2 


40.6- 


- 101.6 


8 


240 


Ha + CII 


1999 Sept.,17-Dec, 2 


239.4 


-315.2 


15 


423 


Ha + CII 


2000 March, 28 - June, 23 


432.6 


-519.4 


10 


275 


Ha + CII 


2000 Sept, 5 -Dec, 8 


592.7 


-687.2 


32 


320 


He i ,15876 


1999 March, 2 - April, 24 


40.6 


-93.5 


4 


240 


He i ,15876 


1999 Sept, 17 - Dec, 2 


239.4 


-315.2 


15 


423 


He i ,15876 


2000 March, 28 - June, 23 


432.6 


-519.4 


10 


275 


He i ,15876 


2000 Sept., 14 - Dec, 8 


602.3 


-687.2 


17 


320 



Finally, the spectra were corrected to the stellar rest frame for a 
radial velocity of -12 km s _1 and normalised to the continuum. 

Ritter Observatory, Toledo, USA - The Ritter spectra were 
obtained with a 1 -m telescope, fiber-fed echelle spectrograph, 
and Wright Instruments CCD camera at the University of 
Toledo over a period of two months (September - November) 
in 2000. The spectral resolving power R 26 000, with the res- 
olution element having a FWHM of about 4 pixels. The spectral 
coverage consists of 9 separate 70- A regions in the yellow and 
red. All the exposures were 1 hour in duration. 

The raw frames were reduced with Ritter Observatory's 
standard reduction script under Sun/IRAF 2.11.3. Removal 
of the telluric lines was done with the IRAF task telluric. 
The template spectra used in this context were artificial rows 
of Gaussians constructed from spectra of telluric standard 
stars taken under various conditions. For each spectrum of 
HD 199478, the template that provided the best telluric cor- 
rection was used. However, for a few spectra none of the tem- 
plates in the library were completely successful at removing 
telluric lines. The spectra were then Doppler corrected to the 
heliocentric rest frame and normalised to the continuum. 



2.2. Consistency check 

An important point of any study which relies on observations 
collected at various observatories, with different instruments 
and equipment, is the mutual consistency among the corre- 
sponding datasets. The ideal way to perform a consistency 



2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under contract with the National 
Science Foundation. 



check is to compare strictly simultaneous data collected from 
different places. 

Fortunately, our sample has three such spectra, that were 
taken at the NAO, TO and RO within 8 hours in the same night 
(Sept. 17, 2000). Using these spectra we checked for possible 
systematic differences in continuum and wavelength calibra- 
tions. The results obtained indicate that the wavelength cali- 
bration of the NAO and the TO spectra agree perfectly while 
the Ritter spectrum shows a one pixel systematic shift to the 
red. 

On the other hand, and as regards photometric calibration, 
the Ritter spectrum (corrected for the shift of one pixel) fits 
quite well (within the noise) the NAO spectrum, while the rela- 
tive fluxes between 6561 to 6584 A in the TO spectrum are up 
to 3% stronger. Since the three spectra are not strictly simulta- 
neous and since line profile variations on a shorter (hours) time 
scale cannot be excluded, it is not currently possible to judge 
to what extent the established differences in the fluxes redward 
of the emission peak of Ho- might be caused by imperfect con- 
tinuum rectification or are due to real variability in the wind. 

Thus, differences of ~9 km s _1 in velocity scale and up to 
3% in relative flux cannot be excluded in our spectral time- 
series (but see next section). 

3. Photospheric variability 

3.1. Photometric evidence 



Recent results (iPercy et al. H 2008) indicate that the photomet- 
ric behaviour of HD 199478 is characterised by continuous 
irregular/multi-periodic variations with an amplitude of about 
0.15 mag on a time-scale of 20 to 50 days. In some obser- 
vational runs colour variations of up to 0.05 mag, in phase 
with the light curve, have been also observed while in others 
no colour variations were detected above the corresponding er- 
ror. In these properties HD 199478 is similar to other OB SGs 
which are known to be photometrically variable and show small 
amplitude microvariations in the visual, with little c olour varia 



tions, on a time scale from days to months (see e.g. lAerts et al 
1999Hvan Genderen 11200 lllMathias et all200lh . 



3.2. Photospheric variability traced by C II doublet and 
He I 6678 lines. 



Following Mar kova & Valchev (2000), we used the absorp- 
tion lines of C n /U6578.03, 6582.85, He i ,16678.15 and 
He i ,15875.67 to probe the deep-seated variability and pho- 
tospheric structure of HD 199478 during the period covered 
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by our observations. To improve the internal consistency of 
the wavelength scale in the extracted spectra from different 
observatories, which is of crucial importance for the purposes 
of the time-series analysis, the C n and He i -16678 line profiles 
were realigned using the diffuse interstellar band at -16613.6 
as a fiducial. Similarly the interstellar line of Na i D -15889.95 
was aligned for our study of profile changes in He i -15876. 
Following these adjustments, we estimate that the velocity 
scale local to each line profile is stable to 1-2 km s -1 . We 
are also confident that the C n lines are not severely affected 
by large fluctuations in the outer red wing of Ha. For the 
photospheric analyses the C n lines were normalised to a local 
continuum assigned (using a low-order polynomial) between 
-1-16570 to 6590A. 



Radial velocities of the selected He i and C n lines were 
measured by fitting Gaussian profiles in the knowledge that 
these lines are generally very symmetric. For the 2000 datasets 
the following estimates of the mean radial velocity and peak-to- 
peak amplitude were derived: -1.9+3.9 km s -1 and 16 km s -1 
for C ii -16578; -2.0 +3.7 km s 1 and 12 km s 1 for He i -16678 
and +2.8+3.7 km s -1 and 13 km s -1 for He i -15876. Variations 
of ~ 15% in the total equivalent widths of the lines were also 
established. There is a tighter correlation between the strength 
and velocity changes seen in C n and He i -16678 than between 
either of these lines and He i -15876. 

The sampling rate of the 2000 (and the 1999) dataset is 
rather uneven and short-time series secured over a few days are 
separated by data gaps of between 1 to 3 months. This makes 
the search for periodic sign als rather more unce rtain. We ap- 



plied the CLEAN method ([Roberts et al.lll987l) to the radial 



velocity measurements (using a gain of 0.5 and 200 iterations). 
The power spectra, where the features of the window function 
have are deconvolved using the discrete Fourier Transform, are 
shown in FigQ] for the 1999 Cn and the 2000 He i -16678 and 
C ii data. 

Clearly, there is no strictly periodic signal present in the 
photospheric lines of HD 199478 that remains coherent be- 
tween 1999 to 2000. There is instead some indication that the 
absorption lines are semi-modulated in their central velocities 
over time-scales of ~ weeks to months. The only signal in the 
2000 power spectrum that is consistent between C n and He i 
-16678 is at a frequency of ~ 0.085 days -1 , i.e. a period of 
~ 11.7 days. Interestingly, in 1999 this modulation is essen- 
tially absent, but the strongest peak in the C n dataset at 0.0478 
days -1 corresponds to precisely twice the 1 1.7 days period. We 
find no evidence for a 11.7 days or 23.4 days modulation in 
He 1-15876. However, note that the 1999 data sampling is more 
fragmented, with a gap of around 140 days. The central veloc- 
ities of the C ii and He i -16678 absorption lines phased on both 
these periods are shown in Fig [2] 

Despite the high signal-to-noise and spectral resolution of 
our data (Sect. 2.2) there is also no evidence for sub-features 
travelling blue-to-red (prograde) in the absorption troughs of 
the lines, that might for example be identified in terms of the 
presence of low-order non-radial pulsations. 




0.1 0.8 0.3 0.4 

Frequency (d _1 ) 



Fig. 1. Power spectrum (in arbitrary units) for the C n (solid 
line) and He 1 16678 photospheric absorption lines. The arrow 
marks the 'stable' peaks at ~ 1 1.7 days and 23.4 days. 
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Fig. 2. The absorption central velocities of C n (filled circles) 
and He i -16678 (open circles) phased on the modulation time- 
scales identified in the 1999 and 2000 datasets. 



3.3. Comparison between spectral and photometric 
variability 

The 2000 differential uvby photometry obtained by SA as well 
as the 2000 UBV data collected by JP, though not strictly si- 
multaneous, cover the same time period as the corresponding 
spectroscopic data. The Fourier and self-correlation analysis 
of these data indicate the presence of a periodic variation of 
18 ±4 (UBV) to 21+4 (uvby) days with an amplitude of about 
0.15 mag. The colour curve of this micro-variation is blue at 
the maxima and red at the minima of the light curve, thus re- 
sembling a Cyg variations in BA SGs. 
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Fig. 3. — The 1998 (left), 1999 (middle)) and 2000 (right) time-series of Ha shown as one-dimensional plots (top) and two- 
dimensional grey scale images (bottom). All spectra have been corrected for V sys =-12 km s _1 . Velocity scale given with respect 
to the rest wavelength of Ha . 



The estimated photometric period is somewhat larger but 
still consistent (within 3<x) with the 11.7 day period variation in 
radial velocity of C n and He i /16678 photospheric lines. This 
finding strongly suggests that the same physical mechanism 
(based in the stellar photosphere) is most likely responsible for 
the two phenomena observed, which might be identified as the 
signatures of pulsations. 

However, note that the interpretation of the photospheric 



variability of HD 199478 reported here and in [Percy et al. 



(2008) is not straightforward in terms of pulsation. On the one 
hand, radial pulsations are not likely since: first, the period is 
not stable between the observing runs carried out in different 
years and second, with only one exception, the estimated peri- 
ods are longer than the radial fundamental pulsational pe riod, 
^rad.fund ~8 days, as derived by lMarkova & Valchevl (12000). On 
the other hand, the irregular character of this variability is quite 
similar to that obse rved in other late B SGs and A-type stars 
(Kau fer et al.l ll997). A possible origin for these variations, at 
least for stars with M* < 40 M , is the action of non-radial 
oscillation modes excited by the opacity mechanism. 

In this respect, we note that: 

i) a period of about 20 days, as derived from photometric and 
spectroscopic data of HD 199478, is fully consistent with 
the value inferred via the period-luminosity relation for B- 
type variables with e xcited g-mode oscillations (Fig. 2 in 
Waelkens et al.lll998l) : 

ii) on the HR diagram with parameters derived with 
FASTWIND, HD 19 9 478 e l egantl y joints the group of B- 
type SGs studied by iBurkil (119781) for which g-mode in- 
stability is suggested to explain their variability (Fig. 3 of 
Waelkens et alj|l998l) . 



Therefore, non-radial g-mode oscillations might explain the 
photospheric variability of HD 199478. But our data lacks ev- 
idence for travelling blue-to-red (prograde) features within the 
absorption troughs of the lines, which normally indicate non- 
radial pulsational behaviour. 

Clearly, very extended time-series datasets are requisite for 
extracting reliable long period signals from the irregular ab- 
sorption line changes which characterise B SGs. These targets 
lend themselves particularly to modest-sized robotic telescopes 
equipped with high-resolution spectrographs. 

4. Wind variability 

4.1. Ha monitoring campaigns in 1999 and 2000 

In Figure [3] the Ha time-series for 1999 and 2000 are shown 
as two-dimensional gray-scale images. Above each of the 
velocity-time frames the corresponding one-dimensional spec- 
tra are plotted to provide a visual assessment of the size of the 
fluctuations at each velocity bin. Gaps between observations, if 
equal or larger than 1 .0 day, are represented by black bands. All 
spectra have been corrected for the systemic velocity, V sys = - 
12 km s _1 . The zero point in velocity corresponds to the rest 
wavelength of Ha . A similar plot, illustrating the Ha t ime- 
series obtained in 1998 (from Mark ova & Valchev 2000), is 
also provided for completeness. 

Figure[3]demonstrates that the Ha profile of HD 199478 is 
strongly variable, exhibiting a large diversity of profile shapes 
and behaviour patterns. In particular, and as also noted by 
Markova & Valchev (2000), in June- July, 1998 as well as dur- 
ing the first two months of 1999, the profile appeared fully 
in emission evolving from a double-peak morphology with a 
blue component that is stronger than the red one, to a single- 
peaked feature centered almost at the rest frame. Some hints 
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about the subsequent development of this feature to the red 
seem a lso to be present. Th r ee suc h cycles have been identi- 
fied by iMarkova & Valchev (2000) (see left and middle pan- 
els of Figure [3): the first - between HJD 2450968-982; the 
second - between HJD 2450998-1009 and the third - be- 
tween HJD 2451 178-189. Another cycle taking place between 
HJD 2451 217-247 can be now be easily recognised thanks to 
the new observations in March 1999. This finding implies that 
the variability pattern described above is relatively stable (over 
at least 9 months) with a characteristic time-scale of about 15 
days and a possible re-appearance after one month or longer. 

A new variability pattern is revealed by the 1999 and + 
2000 observations, where Ho- appears not only in emission, 
but also in partial or complete absorption. In particular, ono 
HJD 2451 293 (April 24, 1999) in addition to the blue-shifted '5 
emission (V r =-75km s ) a slightly red-shifted absorption fea- Q 
ture (V r =+40 km s _1 ) has appeared giving rise to a reversed 
P Cygni profile. The latter persisted for at least 8 days, grow- 
ing slightly stronger in intensity. 

About five months later, namely on HJD 2451439 (Sept. 
17), Ha is seen fully in emission again, though with a weak dip 
at about +60 km s _I , which makes the profile appear double- 
peaked with a blue component being much stronger than the 
red one. This configuration was preserved for at least 7 days, 
i.e. up to HJD 2451 449. On HJD 2451 475 the absorption dip 
is missing but one month later it appears again, stronger than 
before, and persists for at least 6 days (between HJD 245 1 509- 
515) during which time the profile again looks like a reverse 
P Cygni. Interestingly, the second appearance of the dip is ex- 
actly at the same position as the first one suggesting the same 
physical origin for both events. 

The 2000 observations importantly revealed (right panel of 
Figure O the presence of another unusual event during which 
Ha changes suddenly and drastically from pure emission to 
pure absorption and back to pure emission. By chance, the dis- 
tribution of the available observations in time was quite good 
allowing the development of this spectacular event to be fol- 
lowed in more detail. 

In Figure H] one can see that before the onset of the high- 
velocity absorption (HVA) event, Ha appeared fully in emis- 
sion developing from a double-peaked to a single-peaked mor- 
phology and strengthening slightly with time (HJD 245 1 792- 
816). On HJD 2451 827, in addition to the emission a localised 
high-velocity (V r =-150 km s ) absorption extending from -68 
to -250 km s is present making the profile appear P Cygni- 
like. Over the next 9 days the P Cygni feature evolves into dou- 
ble absorption with central emission where the blue component 
is significantly stronger and wider than the red one. Two weeks 
later, (HJD 2451 863) the morphology of the profile is still the 
same though the blue component is weaker and narrower while 
the red one has apparently strengthened becoming somewhat 
wider. Subsequently the two absorptions are fading in parallel 
and disappear completely on HJD 2451 877. Our observations 
further suggest that the absorption event seen in the Ha time- 
series of HD 199 478 between HJD 245 1 827 - 876 may not be 
unique. Indeed, readers should note that on HJD 2451 632 (see 
right panel of Figure [3} Ha has also appeared as a double ab- 
sorption feature. Unfortunately, due to poor temporal coverage 
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Fig. 4. 
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HVA event observed in the Ha data of HD 199478 in 
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Fig. 5. Top: Ha profile of HD 199 478 illustrating the behaviour 
of pure wind contribution to the line (in velocity space). 



the time development of this feature cannot be followed, but 
given the similarity in the morphology of this profile and the 
one taken, e.g., on HJD 2451 853 we are tempted to speculate 
that about 6 months earlier an absorption phenomena similar 
to the one recorded in September - October 2000 may have oc- 
curred in this star. 



N. Markova et al.: Wind structure of late B supergiants. 



7 



4.2. Wind contribution 

To probe further the nature and origin of peculiar emission and 
HVAs in the Ho- data of HD 199478 we normalised the ob- 
served Ha profiles to a constant photospheric profile com- 
puted by means of the F ASTWIND code with parameters from 
Markova & Pulsl d2008h (given also in Table©. 

In Figure [5] the Ha profiles from the 1999 and 2000 data- 
sets are shown (in chronological order from the bottom up- 
wards). Profiles from runs separated by large gaps are grouped 
together where periods without significant Ipv are represented 
by a single averaged profile so as not to confuse the figure. 

The following points are immediately apparent from these 
plots: 

• Outside the HVA events the wind contribution can be as- 
signed to two components: (i) a strong emission feature which 
is either symmetric with respect to the stellar rest frame or 
shows weak blue-to-red asymmetry with a blue wing being 
more extended and stronger than the red one and (ii) localised 
emission bumps with variable position which more likely give 
rise to the established VI R variations. 

• During the HVA episodes (right panel, bottom and middle 
groups) in addition to the extended blue-shifted absorption an 
emission component also exists, i.e. the wind does not only ab- 
so rb but also emits Ha p hotons, contrary to the cases described 
by lKaufer et al. ( 1996bh . where the HVAs are not accompanied 
by unshifted emission. 

• The red-shifted absorption seen occasionally in Ha, is more 
likely of wind origin and suggests the presence of matter infall 
at the base of the wind. 

From the properties outlined above one might conclude that 
the envelope of HD 199478 consists of (i) a spherical com- 
ponent where the physical conditions favour only processes 
which produce emission in Ha and (ii) localised large-scale 
wind structures with matter infall and outflows, where emis- 
sion and/or absorption can originate. 



4.3. Comparisons of high-velocity absorption event in 
HD 199478 to those in other late-B SGs 

Comparison of our Figure H] with similar results from 
Kaufe r et al. dl996albh showed that the spectacular absorption 
event seen in Ha of HD 199478 is qualitatively similar to the 
HVAs observed in HD 34085 (B8 la, j3 Ori), HD 91 619 (B7 
la) and HD 96 919 (B9 la), though with one exception: in our 
data-set the blue and the red-shifted absorption components do 
not merge to form an extended blue-to-red absorption, as is the 
case of the objects Kaufer et al. studied, but instead occur par- 
allel to each other (though we accept the caveat that a more 
intensive and extended dataset is ideally required). 

That the spectacular phenomena of HVAs in Ha have been 
observed so far in 4 late B SGs with peculiar emission in Ha de- 
serves special attention since it might indicate some fundamen- 
tal pr operty of thei r stellar winds. With this in mind we fol- 
lowed |Kaufer_etalJ (11996bl) and measured the main properties 
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Fig. 6. Top: Total equivalent width of the Ha line profile of 
HD 199478 during the 2000 HVA event as a function of time. 
The horizontal lines represent the "undisturbed" mean equiva- 
lent width of Ha (fully-drawn) and its corresponding standard 
deviation (dotted). 

Bottom: Time behaviour of the differential y-magnitudes dur- 
ing the HVA event 

Table 3. Properties of HVAs in Ha observed in 3 late-B SGs. 
Data for HD 34 085 and HD 96 919 are taken fromlKaufer eTaT 



( 1996b). All velocities, measured with respect to the stellar rest 
frame, are given in units of the corresponding wind terminal 
velocities (from Table 0J. 





HD 34 085 


HD 96 919 


HD 199 478 


Signature 


1994 


1995 


2000 


MJD of max. blue depth 


2449493 


2449792 


2451836 


Max depth in % of cont 


20 


70 


49 


Vel. of max. depth 


0.40...0.60 


0.29...0.43. 


0.37...0.76 


Blue-edge velocity 


0.79... 1.20 


0.43...0.64 


0.68...1.39 


Red-edge velocity 


0.34...0.52 


0.32.. .0.48 


0.41...0.84 


Rise time of W,t [d] 


11 


21 


22 


Decay time of W A [d] 


20 


46 


33 


Duration of event [d] 


40 


90 


55 



3 Such an approximation is legitimate since the observed photo- 
spheric variability in HD 199 478 is indeed very weak (see Sec. l3.2t 



(such as, e.g., relative intensity, position and blue- and red-edge 
velocities) of the 2000 HVA in Ha of HD 199 478 at the time 
of its maximum intensity. Based on the time evolution of the 
Ha total equivalent width (measured by integrating the flux be- 
tween 6554 and 6570 A) we determined the total duration and 
the rise and decay times of this event (Figure [6] top panel) 

The estimates thus derived are listed in Table [3] to- 
gether with similar data for HD 34085 and HD 96 919 (from 
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Kaufer et aL 1996b). Note that since in our spectra the blue- 
shifted component of the HVA event never merges with the 
red-shifted one, the red-edge velocity (given in Table|3]in bold) 
does n ot refer to the extend ed blue-to-red absorption (as is the 
case of Kaufer et al . 1996b) but instead corresponds to the red- 
absorption component itself. Note also that due to the large un- 
certainties in the adopted terminal velocities (see next section) 
the normalised velocities sometimes exceed unity. 

Compared to similar events in HD 34085 and HD 96919, 
the HVA seen in HD 199478 is of intermediate duration and 
strength. Its time development is roughly consistent with re- 
sults from Kaufer et al., which show rise times that are smaller 
than the decay times. □ Thus, the duration of a HVA event 
seems to depend on its maximum strength (stronger maximum 
absorption - longer duration), while its development in time 
(rising time vs time of decay) appears to be independent of this 
parameter. In addition, the blue-edge velocity and the velocity 
of maximum depth of a HVA event may anti-correlate with its 
strength, i.e. stronger features tend to reach maximum depth 
at lower velocities, being less extended in velo city space than 



weake r ones. Furthermore, and as also noted bv llsraelian et al 
(119971) . the maximum positive velocity of a HVA is always 



lower than the corresponding maximum negative velocity. (Due 
to the limited number of stars these results can only be regarded 
as suggestive and they have to be confirmed with improved 
statistics.) 

Finally, we note that the photometric behaviour of 
HD 199478 during the 2000 HVA event in Ha provides ten- 
tative evidence that at the onset of the event the star was about 
one magnitude fainter than at the moment of maximum line 
absorption (Figure[6] bottom panel). 

5. Stellar and wind parameters of late-B stars 
which have exhibited HVAs in Ha 

To help understand the nature and o rigin of the Ha varia bility, 
especially the appearance of HVAs, iKaufer et al.1 (Il996ah have 
determined the fundamental parameters of their sample stars 
employing: 

i) the Azzopardi spectral type-absolute magnitude and ob- 
serve d Hy equivalent width-absolute magnitude calibra- 
tions dAzzopardilll98ll) : 

ii) the Schmidt-Kaler spectral type-bolometric correction 
and spectral type -effective temperature calibrations 
dSchmidt-Kalerll 19821) : 



iii) the evolutionary tracks of ISchaller et alJ 41992) to derive 
stellar masses; 

iv) the maximum half-width of the TVS of two isolated ab- 
sorpt ion lines to esti mate projected rotational velocities; see 
also dReidetal.il 19931) 

v) the UV resonance lines Mg 11/1/12795,2803 to determine 
wind terminal velocities, v ro . 

By means of these parameters the authors have subse- 
quently constrained the true rotational periods (by means of 



P 10t /sini i and P rot; break) an d evaluated the radial fundamental 
pulsational periods of the stars. The analysis led the authors to 
suggest that rotation plays an important role in determining the 
properties of the Ha variability including the HVA events. 

Prior to the use of currently available state-of-art model at- 
mosphere codes, the approach used by IKaufer et al.l (11996a) 
(with its well-known weaknesses and uncertainties) was the 
only one to permit the basic parameters of hot stars to be deter- 
mined. 

The situation has changed drastically since then and stel- 
lar and wind parameters of hot stars can now be derived 
with relatively high precision using the methods of the quan- 
titative spectral analysis. The outcomes of such analyses, 
performed by means of the present day N LTE, line bla n- 



keted model atmosphere codes (e.g. CMF GEN lHillier & M iller 
19981 and FASTWIND IPuls et al. Il2005l) have unambiguously 
showed that the newly derived stellar and wind parameters 
can significantly dev iate from their earlier determinations (e.g., 
Martins et al. 20051 and references therein for O stars and 
Crowther et alj|2006t iMarkova & Puldl2008llSearle et~ai]|2008 
for B stars) 

With this in mind and given the limited number of late B 
SGs with reliably deter mined stellar and wind parameters (see 
Markova & Pulsll2008l and references therein), we decided to 
re-determine the parameters of the Kaufer et al. late-B SGs 
with HVAs in Ha , using optical spectra kindly provided by 
Otmar Stahl and employing the one of the latest version of the 
FASTWIND code. This way a homogeneous data base for late 
B SGs sharing similar empirical properties in Ha would be cre- 
ated, which might be easily extended in the future. 

To perfor m our analysis w e follow ed the strategy outlined 
in detail by IMarkova & Pulsl ( 2008 ). In particular, effective 
temperatures, T e g , were estimated from the Silicon ionization 
balance, fitting the Si II doublet at 4130 A and the Si III triplet 
at 4552 A and adopting a solar Silicon abundan ce (log (Si/H)= 
-4.45 by numbeid, cf. Grevesse & Sauval fl998 (and references 
therein), and a microturbulent velocity, v m i c , approp riate for 
the corresponding spectral type (IMarkova & Pul|2008). Since 
in all objects the blue wing of Hy seems to be affected by blue- 
shifted wind emission (similar to the one seen in Ha ) surface 
gravities, log g , were derived fitting the wings of H6 . The ac- 
curacy of these estimates is ±500 K in T e g and ±0. 1 in log g . 



Projected rotational and macroturbulent velocities, vsin; 
and v mac , - these parameters were determined from Mg II 
A 4481 employing the Fo urier technique developed by 
Simon-Diaz & Herrerol (12007b . Since this method provides 
only rough estimates of v mac the latter have been addition- 
ally adjusted during the fitting procedure (if nece ssary) to im- 
prove th e quality of the fits. Radial velocities from lKaufer et alJ 
(1996a) lave been adopted. 



4 This result has to be considered with caution since the exact time 
of maximum depth absorption in HD 199 478 is not known with con- 
fidence due to limited time-series coverage. 



According to latest results dAsplund et aT1 l2005), the actual solar 
value is slightly lower, log (Si/H) = -4.49, but such a small difference 
has no effect on the quality of the line-profile fits. 
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Table 4. Stellar and wind parameters of HD 34085 (J3 Ori) and HP 96919 a s derive d in the present study employing the 
FASTWIN D code. Numbers in brackets refer to similar data from iKaufer et alj dl996al) . Estimates for HD 199478 are taken 
from Markova & Puis (2008). Effective temperatures are given in kK, velocities in km s 1 , time periods in days. Estimates for 
M and logD mom are lower limits. 



star 



sp 



distance Mv 



log L/L Q 



M, /Mo log M 



logD u 



HD91619 


B7Iae 


2.51 


-7.00(-7.99) 


13.9(12.2) 


1.85(1.75) 


63(114) 


5.13(5.42) 


11(27) 


-6.92±0.32 


27.18±0.33 


HD 199478 


B8Iae 


1.84 


-7.00 


13.0 


1.70 


68 


5.08 


9 


-6.73.. .-6.18 


27.33. ..27.88 


HD34085 


B8Iae 


0.50 


-8.31(-7.77) 


12.5(11.2) 


1.70(1.67) 


129(116) 


5.56(5.28) 


31(23) 


-6 47+ - 27 

u '^'-0.15 


27 75+ - 27 

-0.20 






0.24° 


-6.70 






61 


4.92 


7 


-6 96 +0 37 


27 1 +037 


HD 96 919 


B9Iae 




-7.0*(-7.97) 


11.0(10.3) 


1.50(1.50) 


71(141) 


4.82(5.30) 


6(23) 


-7 10 +025 

'■ L -0.13 


27 03 +0 - 27 


star 


v sys 




"mic 


v sin i 


Vbrcak 


P to t, break 


Ptat/sini 


Vesc 


Voo 




HD91619 


-6 


35 


8 


35(60) 


153(212) 


21(27) 


91(96) 


220(261) 


170...330 




HD 199 478 


-12 


40 


8 


41 


134 


25 


84 


191 


170...350 




HD34085 


+ 18 


35 


8 


30(55) 


184(195) 
129 


35(30) 
24 


218(107) 
103 


262(244) 
183 


230...350 
230...350 




HD 96 919 


-24 


25 


7 


30(60) 


112(176) 


32(40) 


120(119) 


160(220) 


250.. .370 





' - HIPPARCOS distance estimate 

- absolute magnitude from the calibration of Humphreys & McElrovl J 19841) 



Stellar radii, R* - These parameters were determined from 
the derived effective t emperatures and de-reddened absolute 
magnitudes (see, e.g.. lKudritzkilfl980l) . The latter were cal- 
culated using standard extinction law with R = 3.1 combined 
with: (i) visual magnitudes, V, and B — V colours from the 
HIPPARCOS Main Catalogu e (I /239/ hip mni „); (ii) intrinsi c 
colours (B - V) = -0.03 from dl-it/patrick & Garmanvl 1 99(j). 
and (iii) distances collected from various sources in the litera- 
ture. 

In particular, for HD 91 619, a member of Car OB 1 associ- 
ation, a distance of 2.51 kpc as provided by Humphr eys! ( 1978 ) 
was adopted. For HD 34 085 (J3 Ori) the situation is a bit more 
complicated. As a mem ber of Ori OB 1 thi s star should be sit- 
uated at about 0.5 kpc ( Humphreys! 1978). Its possible m em- 
bership of the t Ori Rl complex dHoffleit & Jaschekll 19821) re- 
duces the distance to about 0.36 kpc, while the HIPPARCOS 
distance estimate is 0.24 kpc. @ Thus, for B Ori we provide 
two entries as upper and lower limits to the distance to account 
for all possibilities. For HD 96919, which does not belong to 
any cluster or asso ciation, an absolute magnitud e according to 
the calibration of iHumphrevs & McElrov 0984) was adopted. 
We quote typical uncertainties of +500 K in our estimates 
and of +0.4 in My (for members of associations) to 0.5 ma g 
(for stars with My from calibration) (iMarkova & Pulsll2008l) . 
The error in the stellar radius is dominated by uncertainties in 
M v and is of the order of Alog = +0.08. ..0.10, i.e., less than 
26% in R+ . 



Luminosities, log L/L Q , and stellar masses, M* , - these 
values were determined from the corresponding T^and 
R+ values and the "true" surface gravities, respectively. LI The 



typical uncertainties of these estimates are Alog L/L G - 
to 0.22 and AlogM* =+0.19 to 0.22. 



t0.17 



6 Below 0.5 kps the HIPPARCOS distance estimates are generally 
accepted as reliable. 

7 "True" gravity results from the observed gravity corrected for the 
centrifugal acceleration (=(v sin i ) 2 )/ R* ). Due to the lower v sin i of 



Terminal wind velocities^ , - values for B Ori and 
HD 96919 have been determined from the blu e-edge of the 



Mg II resonance lines at AA 2795, 2803 by Kaufer et al 
dl996al) . However, the authors note that due to the absence of 
sharp blue edges, their estimates have to be considered as lower 
limits only. Thus, in these two cases we adopted the Kaufer 
et al. estimates of v^, (i.e. -230 km s _1 for HD 34085 and - 
250 km s _1 for HD 96 919) but assumed an asymmetric error of 
+50% to allow for a rather large uncertainty towards higher val- 
ues. For HD 91 619, since no Vqq estimate w as found in the lit- 
erature, we followed IMarkova & Pulsl (12008b and adopted v ra = 
v esc = 220 km s _1 assuming an asymmetric error of -25/+50%. 

Mass-loss rates, M , and velocity exponent ft - In the case 
of strong undisturbed winds (Ha in emission) M and B can be 
estimated, with relatively high precision, from the best fit to 
the red wing and the peak emission of the Ha profile, respec- 
tively. Indeed, Ha is in emission in the spectra of HD 91619 
and HD 96919 available to us. However, and as also shown in 
Figure|7] such profiles cannot be reproduced in terms of spher- 
ically symmetric smooth wind models since at this temperature 
regime the models predict profiles in absorption partly filled 
in by wind emission. Due to this reason only lower limits to 
Mare derived for our sample stars, with f3 ranging from 0.8 
to 1.5, except for HP 96919 where an upper limit to f3 of 1.3 
was adopted (see IMarkova & Pulsll2008l) . Under these circum- 
stances any excess emission seen in Ha should be attributed 
either to deviations from spherically symmetric wind approxi- 
mations and/or to processes different from recombination (see 
sect. [6l). 



the sample stars this correction is generally small, between 0.01 to 
0.03. 
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Fig. 7. Examples of typical Ho- profiles observed in the four sample stars which have been recognised to show HVA events. 
Over-plotted are synthetic profiles accounting for the contribution of photospheric absorption plus spherically symmetric smooth 
wind emission (for more info see text.) 



The errors in our M estimates (actually in log0 accumu- 
lated from uncertainties in /3 and in R+ are typically less than 
0.30 dex. Having , M and Vco thus determined we finally cal- 
culated the modified wind momentum, £> m om=2voo R* 2 , with a 
typical error AlogD mom between 0.13. ..0.35 dex. 

5.1. Comparisons to parameters derived in other 
investigations. 

Stellar properties derived i n our analysis, togethe r with the 
data for HD 199478 (fro m iMarkova & Pulsl 120081) . are listed 
in Table [4] Estimates from lKaufer et alU 1996ah are also pro- 
vided (numbers in brackets) for comparison. 
Compared to similar data from Kaufer et al.: 

• Our v sin i estimates are about a factor of 0.5 lower. This result 
is easy to interpret since our estimates account for the effects 
of macroturbulence. Contrary to what might be expected, the 
lower values of v sin i did not result in significantly different 
values of the rotational period (see Table[4]l. With our approach 
the decrease in v sin i is almost completely compensated by an 
increase in (see below). 

• Our r e jf and log g estimates are surprisingly similar to those 
in Kaufer et al. This finding indicates that in these particular 
cases the effects of line blanketi ng are not too large (as also 
shown bv lMarkova & Pulsl 120081) and that in the Kaufer et al. 
approach the effects of the larger stellar masses have been to a 
large extent compensated by the larger radii. 

• Our R+ and log L/L G are significantly lower than those of 
Kaufer et al. The reason for this is the large difference in the 
adopted absolute magnitudes: according to our estimates the 
three stars are about one magnitude fainter than estimated via 
the Azzopardi calibrations (apart from HD 34085 (J3 Ori), first 
entry, where there is an alternative possibility of My being half 
a magnitude brighter). Such large differences in My cannot be 
solely explained by differences in the methods used, but might 
indic ate instead problems w ith the Azzopardi calibration (see 
also lMarkova & PulslEooi) . 

• In general, our estimates of M* tend to be lower (by about 
0.39 to 0.58 dex) than those estimated by Kaufer et al. via 
the evolutionary tracks (apart from the first entry for /3 Ori, 
where a value of 0.12 dex larger was derived). This finding 
is consistent with similar results from previous investigations 



(Cro wfher et alj 120061: iTrundle et al l 120051: IMarkova & Puis 
2008) which also indicate a mild "mass discrepancy" for B 
SGs. Indeed, the differences established here are somewhat 
larger, but in their study Kaufer et al. used evolutionary tracks 
whi ch do not take into account the effects of stellar rotation, 
e.g.. lSchalleret~ai1(ll 992). while in all other studies evolution- 
ary tracks from lMeynet & Maeder ( 2000l) (with stellar rotation) 
have been used instead. 

On the other hand, the stellar radii and luminosities de- 
rived injhe_rjjreseji^tudyare consistent with similar results 
from IMarkova & Pulsl (120081) . In particular, they fit very well 
(within the existing scatter) the log Teff-7?* and log Teff- 
logL/L distributions derived by these authors for Galactic B 
SGs. The only real outlier is HD 34085 (J3 Ori), first entry, 
where larger deviations (~ 0.6 dex in log L/L Q and a factor of 
2 in Ri, ) have been established. Since the log g -value of this 
star is typical for a l ate SGs we sugge st that the photometric 
distance provided by lHumphreysl (119781) is overestimated. 

Finally, strong agreement (within the corresponding er- 
rors) was found between our T e g , log g and v m j c estimates 



of HP 34085 (J3 Ori) and those derived by iPrzvbilla et al 
d2006l) via a hybrid non-LTE technique (r eff =12 000+200 KK, 
logg =1.7+0.1, v m j c =7±l km s" 1 ) Good agre e ment was 



also found with the results from llsraelian et al. 
HD 34085 (r eff 



19971) for 

13000, log g =1.6 and v mic =7) obtained 
via the NLTE unb lanketed plane-parallel hydrostatic code 
TLUSTY dHubenvll988l) . The latter indicate that at these tem- 
peratures the effects of line blocking /blanketing are small as 
also found bv lMarkova & Pu is (2008) and that in the particular 
case of this star the wind effects also seem to be minimal 0. 



6. Discussion and conclusions 

Extensive monitoring campaigns of severa l late-B SGs, namely 



HP 199478 (present study as well as Markova & Valchev 
2QOCj) an d HP 91 619, HD 34 085 an d HP 96919 dKaufer et al 



199 6aR[l997t llsraelian et alll997h . indicate that their Ha pro- 
files exhibit quite similar peculiarities, consisting of a double- 
peaked emission with V/R variation and occasional episodes of 
strong absorption indicating simultaneous mass infall and out- 



8 (Although the somewhat lower value of log g of Israelian et al. 
more likely results from the neglect of the wind effects) 
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flows. Such line signatures cannot be reproduced by conven- 
tional (i.e. non-rotating, spherically symmetric, smooth) wind 
models, which instead predict profiles in absorption partly 
filled in by emission for SGs at this temperature regime. 

Discrepancies between observed and predicted Ha profiles 
have been also established for many O and early B SGs, where 
this finding was usually interpreted as an indication for devi- 
ations from the adopt ed spherically symmetric, smooth wind 
appro ximations (e.g.. iMorel et al. 2004 : Markova et al. 2004, 
2005 and references therein). 

Following this reasoning, axially symmetric envelopes, 
modulated, at least in the inner parts, by co-rotating weak 
magnetic structures have been assumed to explain the appear- 
ance and kinematical properties (e.g. double-peaked morphol- 
ogy with V/R variations) of the peculiar Ha emission in the 
spectra of the four late-B S Gs noted above dKaufer et al.ll9 96a; 
Mar kova & Valchevll2000l) . 



In addition, to account for the sudden appearance of HVAs 
in Ha and their development in time (e.g. the fast rise over a 
large velocity range, the lack of unshifted line emission, and the 
mooted re-appearance over a rotational time-scale), episodic 
and azimuthally extended, density enhancement s in the form 



of co-r otating spirals rooted in the photosphere dKaufer et al 



996bh or clos ed magnetic loops similar to those in our Sun 



dlsraelian et al . 1997) were also suggested 

It is generally mooted that non-radial pulsations (NRPs) 
and surface magnetic spots may equally be responsible 
for creating large-sc ale inhomogeneities in hot star winds 
dFullerton et aDl l996>. However, despite some progress (e.g., 



Kaufer et al. 2006) no convincing evidence of a direct relation 



between the time-scale of a given cyclical (wind) Ipv and the 
predicted time-scale of recurrent surface features du e to a spe- 
cific pulsation mode, has been derived to date (see iTownsend 



2007 and references therein). 



Note in the particular case of the four stars discussed here 
non-radial pulsations due to g-modes oscillation s were sug- 



gested to ex plain absorption Ipv in their spectra ( Kau fer et al 



19971 iMarkova & Pulsl l2008). This possibility is partially sup 



ported by the present results, which indicate that on the HR dia- 
gram, and for parameters derived with FASTWIND, these stars 
fall exactly in the region occupied by known variable B SGs, 
for which g-modes instability was suggested. Also, the photo- 
metric variability of HD 199478 seems to be co nsistent with 
a possible origin in terms of g-mode oscillations ( Per cy et al. 



a pos 
2008 



2008). Thus, it seems very likely that the four late-B SGs in our 
sample are non-radial pulsators. Although no clear evidence of 
any causality between photospheric and wind (as traced by Ho-) 
variability has been seen so fa r for these objects (present study 
as well as iKaufer et al.lll997l) . one might speculate that their 
winds are perturbed due to pulsational instability, with specific 
signatures seen in the behaviour of Ho-. 

An alternate possibility is that magnetic fields could 
be responsible for the appearance of large-scale structures 
and wind asymmetries in hot stars. In particular, magneto- 
hydrodynamical (MHD) simulations for stars with moderately 
strong rotation, and for stellar and wind parameters typical 
for O and early B SGs (plus a magnetic dipole aligned to 
the stellar rotation) showed that depending on the magnetic 



Table 5. Magnetic field strength, B (in G), required to get 
an equatorial confinement with simultaneous mass infall and 
outflows around each of our targets. The Keplerian, R&, the 
Alfven, Sa> and the escape, Re, radii ( in units of R+ above 
the ph otosphere) are calculated following lOwocki & ud-Doula 
(2003). The HD 34085 estimates correspond to the 2nd entry 
in Table H 



Parameters 



91619 



199 478 



34085 96919 



Ra 
Re 

B (R K <R A < R E ) 



0.37 
0.83..3.73 
3.77 
5.. 100 



0.45 
0.77..4.59 
4.58 
5..180 



0.38 
0.83..3.79 
3.81 
5..105 



0.42 
1.0..4.13 
4.18 
5..85 



spin-up, an equatorial compression, dominated by radial infall 
and/or outflows can be created, with no apparent tendency 



to form a steady, Keplerian d isk ( Owo cki & ud-Doul a 2003 



ud-Doula. Owocki & Townsendll2008h 

Indeed, due to the lack of strong convection zones asso- 
ciated with hydrogen recombination, normal (i.e. without any 
chemical peculiarities) hot stars are generally not thought to be 
magnetically active. However , theoretical considerations (e.g. 
Cassinelli & Macgregot 2000l) supported by more recent ob 



serva tions dBvchkov. Bvchkova and Mad ei 2003; Hubrig et al 



120071) indicate that this may not necessarily be true and that 
relatively strong, stable, large-scale dipole magnetic fields are 
present in different groups of B stars (e.g. SPB, Be, /3 Cep itself 
etc.) 

Thus, it seems likely that in at least some hot stars mag- 
netic fields can be an alternative source of wind perturbations 
and asymmetries. The potential role of magnetic fields in B 
SGs remains intriguing, especially because it might provide a 
clue to understand the puzzling problem of the simultaneous 
presence of red- and blue-shifted absorption in Ha profiles of 
the four late-B SGs discussed here. 

Guided by t hese perspectives, we emplo yed the scaling re- 
lations given in lOwocki & ud-Do ula (2003) and calculated the 
Alfven, Ra, the Keplerian, R&, and the "escape", Re, radii of 
our targets, using data from Table [4] and fixing the magnetic 
field strength at the values required to create an equatorial con- 
finement. Interestingly, the results listed in Table |3J show that 
in all four cases a very weak dipole magnetic field can effec- 
tively channel the wind outflows, leading to the formation of 
an equatorial compression with simultaneous radial mass infall 
and outflow. 

With this in mind, new MHD simulations for the case of 
mid/late-B SGs have been recently initiated. The preliminary 
results (private communication, Asif ud-Doula) indicate that 
a pure dipole magnetic field of only a few tens of Gauss is 
indeed required to obtain a cool equatorial compression (with 
mass infall and outflow) around a rotating star with stellar and 
wind parameters as derived with FASTWIND for HD 199 478. 
(More detailed information about the outcomes of this study 
will be provided in a forthcoming paper.) 

An obvious advantage of the model described above is that 
it has the potential to at least qualitatively account for some of 
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the puzzling properties of the Ha line of our targets. In particu- 
lar, the sudden appearance of red and blue-shifted absorptions 
might be explained if one assumes, that due to some reason 
the plasma in the infalling or outflowing zones of the compres- 
sion or in both of them (during the HVA episodes) can become 
optically thick in the Lyman continuum and L„, thus forcing 
Ha to behave as a resonance line, i.e. to absorb and emit line 
photons. The kinematic properties of the resulting absorption 
features are difficult to predict from simple qualitative con- 
siderations but it is clear that these properties cannot be dom- 
inated by stellar rotation, but instead will be controlled by the 
physical conditions inside the compression. 

Concerning the interpretation of the peculiar Ha emission, 
the situation is more complicated since such emission can orig- 
inate from different parts of the envelope under quite different 
physical conditions. For example, one can expect that the cool, 
less dense plasma outside the compression will only emit line 
photons (via recombinations), producing pure emission fea- 
ture^) in Ha . 

Also, the cool equatorial compression might contribute to 
the Ha emission, providing the plasma inside the compression 
can occasionally become optically thin in this line. However, 
note that even a plasma that is optically thick in La and Lyman 
continuum can, under specific conditions, produce pure emis- 
sion profiles in Ha (e.g. if collisions dominate the Ha for- 
mation or if due to so me reasons the 2nd an d 3rd levels of 
Hydrogen go into LTE ( Petrenz and Puislll996h l. 

Therefore, very weak dipole magnetic fields might be re- 
sponsible for creating wind structures in the envelopes of late- 
B SGs. The models derived via MHD simulations seem to 
have the potential to account, at least qualitatively, for some 
of the peculiar characteristics of Ha in the spectra of our tar- 
gets. However, more detailed quantitative analysis is required 
to investigate this possibility further. 

New high signal-to-noise observations to prove/disprove 
the presence of weak magnetic fields can help to clarify the pic- 
ture. Of course, due to the low strength of the magnetic fields 
required one cannot expect to detect these fields directly but 
indirect evidence such as e.g., the detection of X-ray emission, 
abundance anomalies, specific periodic variations in UV res- 
onance lines, interferom etric observations (for more informa- 
tion see Henrichs might also be considered. Note that a 
weak longitudinal magnetic field of about 130+20 cM . might 
have been detected in HD 34085 (J3 Ori) dSevernylll970l) . but 
confirmation is lacking. 

Finally, at the cooler edge of the B-star temperature regime 
pure emission profiles in Ha can be accounted for if one as- 
sumes the winds are clumped. Indeed, a spherically symmet- 
ric, clumped wind will mimic wind densities higher than the 
actual ones, thus giving rise to strong line emission, similar to 
that in O stars. Such winds may also give rise to wind absorp- 
tion, providing some of the clumps are optically thick in Ha. 



9 (Depending on the size of the Ha forming region emission may 
not appear in the spectrum.) 

10 Given that only one spectral line has been used for these mea- 
surements an error of 20 Gauss seems somewhat unrealistic and likely 
represents only a lower limit. 



Detailed numerical simulations and line formation calculations 
are required to discriminate between the different possibilities. 
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